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Abstract. - In semiclassical studies of systems with mixed phase space, the neighbourhood 
of bifurcations of co-dimension two is felt strongly even though such bifurcations are ungeneric 
in classical mechanics. We discuss a scenario which reveals this fact and derive the correct 
semiclassical contribution of the bifurcating orbits to the trace of the unitary time evolution 
operator. That contribution has a certain collective character rather than being additive in 
the individual periodic orbits involved. The relevance of our observation is demonstrated by a 
numerical study of the kicked top; the collective contribution derived is found to considerably 
improve the semiclassical approximation of the trace. 



Periodic-orbit theory |jj aims at the semiclassical evaluation of energy levels and relates the 
trace of the unitary time evolution operator to periodic orbits of the corresponding classical 
systems. A recent semiclassical study Q based on periodic-orbit theory was devoted to the 
neighbourhood, in the space of control parameters, of classical bifurcations. A collective 
treatment of the bifurcating orbits was found necessary, and even more the inclusion of 
predecessors of such orbits which live in complexified phase space and were termed ghosts. The 
existing semiclassical studies ^, |j, ]5| focus on the generic bifurcations in the classification 
of Meyer || . The purpose of this Letter is to demonstrate that even the neighbourhood of a 
bifurcation which is ungeneric to the classical system is felt semiclassically and necessitates a 
collective treatment. This implies that collective contributions of this kind will constitute a 
basic ingredient in a semiclassical trace formula for systems with mixed phase space. 

The situation we have in mind manifests itself in certain sequences of bifurcations. Sadovskh, 
Shaw, and Delos @] found that such sequences can be explained by normal-form theory || and 
mentioned the importance of their observation for semiclassical studies of systems with mixed 
phase space. We go one step further and derive the explicit collective semiclassical amplitude 
of a group of orbits involved in the bifurcations. 

We shall concentrate on a certain sequence of bifurcations, that of a period tripling and the 
tangent bifurcation of the satellite. In the diamagnetic Kepler problem considered in this 
scenario was not found. 

The semiclassical expression to be established will be tested for the periodically kicked top 
M . A period-to-period stroboscopic description of the quantum mechanics involves the unitary 
Floquet operator 

k k 
F = exp(-i ^ Z +l Jl - ip z J z )exp(-ipyJ y ) exp(-i^-^-y - ip x J x ) . 

We here encounter angular-momentum operators J x ,y,z obeying the commutator relation 
[Jk, Ji] — itijkJk- The Hilbert space dimension is fixed as 2j + 1 (with j the good quantum 
number from J 2 = j(j + 1)), and the semiclassical limit is reached by sending j — > oo. The 
classical phase space is the sphere 3 2 + 1)) = 1. The parameters pi may be interpreted as 
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(a) k=3.524 (b) k=3.53 (c) k=3.545 (d) k=3.57 



Fig. 1. - Section of the phase space displaying the sequence of bifurcations. In (a) a stable orbit is 
surrounded by its stability island. From (a) to (b) two satellites appear via a tangent bifurcation, 
the inner one approaching the central orbit in the period tripling (c) and finally re-emerging on the 
opposite side (d). 

angles of rotation while the ki characterize nonlinear rotations sometimes called torsions. To 
steer the system towards bifurcations, we hold the pi at fixed values and vary k = k z = 10k x ; 
for k = the system is integrable while for k = 5 it displays well developed chaos. Neither 
time reversal nor any geometrical symmetry is present. 

The situation for the period tripling is visualized in the sequence of phase-space portraits 
in fig. |l|. Initially (a), a stable periodic orbit of period one is surrounded by its stability island. 
At a certain value of the control parameter (fc( tan ' = 3.525) two satellites of triple period 
come into existence via a tangent bifurcation. Then the inner (unstable) satellite approaches 
the central orbit (b), collides with it in the period tripling at fc( trip ) = 3.544 (c), and finally 
re-emerges on the other side (d). Astonishing closeness of the bifurcations in parameter space 
seems not exceptional even in the regime of well developed chaos: Starting with k = we found 
the first four period triplings of orbits with primitive period one all accompanied by a tangent 
bifurcation at the pairs (fc< tan ) = 1.853, fc( trip ) = 1.859), (fc( tan ) = 1.971, fc( trip > = 1.975), the 
one mentioned above, and (fc( tan ) = 9.334, fc( trip ) = 9.370). It has to be stressed that no second 
parameter was tuned in the present study to achieve such closeness. However, if one varied 
a second control parameter, the system could be tuned such that both bifurcations happen 
at the same point in parameter space. This means that one deals with the unfolding of a 
bifurcation of co-dimension two. 

The classical stroboscopic map (q,p) — > (q',p') (for one degree of freedom) has the classical 
action S(q',p) as the generating function, with S q < — p' , S p — q (here and in the following 
we denote partial derivatives by indices). Close to the bifurcations (k sa ^.(tan.tnp)^ an j j n 
the neighbourhood of the bifurcating orbits, the action can be approximated by normal forms 
[§). They contain the relevant information about the configuration in phase space and the 
scaling properties of the classical quantities. For the period tripling (without the tangent 
bifurcation of the satellites) one can use action-angle coordinates /, ip with q — \p2J cos ip and 
p = \f2I sin ip and has 

S(I,ip') = If' + el + al 3/2 cosSip' . (1) 

Here e = k — fc( trip ) is the control parameter; at e = the period tripling takes place. 
The coordinates of the satellite periodic orbits are given by Si = if', = I; hence, the 
(/^-coordinate obeys — 3a/ 3 / 2 sin 3ip' = 0. Since a threefold symmetry is implied by this equation 
it suffices to consider the second equation Sj = on the g-axis after switching back to the 
coordinates p, q, yielding eg + ^7g a< 7 2 = 0. This describes the central orbit at I = and the 
unstable satellite at distance I = g 2 /2 = 4e 2 /(9a 2 ). 




H. SCHOMERUS: SEMICLASSICAL INTERFERENCE OF BIFURCATIONS 



3 



The normal form originally stems from a Taylor-Fourier expansion of the full action around 
the central orbit, followed by a certain rectification procedure. While this series generally does 
not converge one still can attempt to include higher order terms to explain certain phenomena 
(see also [|), as was argued for the inclusion of further satellites by the authors of fit]. 

The next-order term is bl 2 , and the extended normal form is 

S(I, p') = lip' + el + al 3/2 cos Zip' + bl 2 . (2) 

The (^-coordinate of the periodic orbits once more obeys — 3a/ 3 / 2 sin3ip' = 0. On the g-axis 
they now satisfy eg + "^g a ° 2 + bq 3 = 0. This equation has three solutions, 



3 a /9 a 2 e 

90 = ' q± = -w = 2-b ± VT2¥-b- 

One in fact sees that the inclusion of the next-order term implies the existence of a further 
satellite. We will denote the central orbit by the index and the satellites by ± as it is done 
here for the coordinates. At g( tan ) = 9a 2 /(326) the satellites undergo a tangent bifurcation 
and for e/b > 9a 2 /(326 2 ) both satellites are ghosts (cf. Fig. g (a)). For < e/b < 9a 2 /(326 2 ) 
both satellites are on the same side of the central orbit (Fig. [I] (b)), while after the period 
tripling (e = 0) they lie opposite to each other, Fig. |l| (d). The situation approaches that of 
a broken torus as e/b — * — oo. When a second parameter is varied to achieve a = e = both 
satellites are contracted onto the central orbit in a co-dimension two bifurcation. 
Our semiclassics for maps starts with the expression || |l0|, [Ilj 

^"=//^i«r-p{^<"'-«-^} (3) 

of the trace of the Floquet operator F as an integral over phase space. It involves the classical 
action S^ n '(q',p) generating the n-th iteration of the classical map and a topological phase 
expressed by the Morse index v. The orbits give a semiclassical contribution C"( clustCT ) to tr F n 
if n is a multiple of their primitive periods , related to each other by n+ = = 3n °^ . 

The asymptotic behaviour for large e/Ti is found with the stationary-phase approximation 
(spa) 0,0, 

k={o,±y 1 ' 

a sum of three individual contributions, one for each orbit. The 2x2 stability matrix M 
describes the linearized map at the locus of the orbit in phase space. The Maslov indices are 
fj.0 = v — signe, \i a = v, /i_ CT = v — sign b with a — signafc. 

Introducing in (j^) the normal form (^) into the exponent with y/Si ilfi > ~ 1 one obtains 

the collective semiclassical contribution C^cai' ^ of the bifurcating orbits to tr F n in the local 
approximation valid close to the bifurcations (e/Ti small). The coefficients in the normal form 
have to be expressed by the actions Sk of the orbits in order to obtain a contribution that is 
invariant under canonical transformations |j| . 

One can find a uniform approximation which interpolates between the local and the 
stationary-phase approximation by expanding 'Si, v > wl + Bl + CI 2 where the corrections 
come from the higher order terms of the normal form. The uniform approximation found is 

=JJ %^(l + BI + CI 2 )e. P ^(eI + aiy 2 co S 3p' + bI 2 )- i ^ (5) 
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Fig. 2. - The action spectrum |T , ' 3 ^(S , )| 2 as the systems is steered across bifurcations. The peaks at 
S ~ 1 belong to the three orbits which take part in the bifurcations at fc( tan > trlp ) « 1.85, the other 
peak to another cluster of orbits bifurcating at fc( tan > trlp ) m 1.97. 



and is valid for arbitrary combinations of e and Ti as these parameters approach zero. The 
coefficients B and C have to be determined such that (||) is recovered as Ti — > 0. 

In the limit e/6 — > — oo the satellites form a broken torus, well separated from the central 
orbit. To discuss this case it is useful to cast the prefactor into the alternative form 1 + B'l + 
(j/j3/2 s j n g^/ p ar tial integration. One then recovers the uniform approximation found by 
Tomsovic, Grindberg, and Ullmo |l2] ]. 

Let us sketch how one arrives at a form of the integral (|5|) which makes a numerical 
evaluation tractable. The integration over ip' gives a Bessel function, and after proper rescaling 
of I one arrives at an integral of type K(e, b) = Jq(I 3 ^ 2 ) exp[i(el + bl 2 )] dl for the local 
approximation. Expanding the exponential function and with the help of identities for the 
T-function at double and triple argument we obtain a useful representation in terms of a 
double sum, 

I— OO OO / . \ jy 

*b» = yEE { ^{-ib)-^{-Aib)-^ (6) 

n— ??i—0 

/ r (ffi + 3m) (4&/«)" 3/2 r(f + 2 + 3m)\ 
X ^(2m)!r(i + m) (2m + 1)!T (f + m) J 

The two additional integrals in the uniform approximation can be found by taking the deriva- 
tive with respect to ie. 

To investigate the quality of the various semiclassical approximations we now turn to the 
study of the action spectrum, that is the Fourier expansion of tiF n with respect to j = 1/Ti. 
Specifically, we study the function |l4j 

64 

T^(S) = _£y«Str F n (j). 

3=1 

The squared modulus of this function shows peaks at values of S corresponding to values 
of the classical actions of periodic orbits, as can be seen using the asymptotic behaviour (Q) 
of the contribution C*( clustor ) to trF n . Fig. illustrates this for the quantum mechanical 
exact |T^ 3 ^(5)| 2 and values of the control parameter close to the period tripling at fc( trip ) = 
1.859. The peak at S w 1 corresponds to the bifurcating orbits. If one steers away from 
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Fig. 3. - Quantum mechanically exact and semiclassically approximated peak heights in the action 
spectrum at the values of the classical actions of the three orbits which bifurcate at k « 1.85, plotted 
as a function of the control parameter k. Both the uniform and the local approximation work well 
close to the bifurcations; the uniform approximation remains valid even far away from the bifurcations. 
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Fig. 4. - Peak height of the central orbit as in fig. |3|. Both the co-dimension-one approximation as well 
as the approximation that treats all orbits as isolated only work well far away from the bifurcations. 



the bifurcations the peak is resolved into three distinct ones at the classical actions of the 
bifurcating orbits. It is instructive to calculate the height of the peaks both quantum- 
mechanically exact as well as semiclassically as k is varied across the bifurcations. Fig. || 
displays the result obtained with the uniform and the local approximation. Close to and on 
the left of the bifurcations both collective contributions are in such good coincidence with the 
exact result that the difference is hard to resolve. Far away from the bifurcations the local 
approximation starts to fail while the uniform approximation remains valid. 

Fig. |^ reveals that the three orbits necessarily have to be treated together close to the 
bifurcations. For k > fc( tnp ) = 1.859, the co-dimension-one approximation treats the outer, 
stable satellite as isolated and the central orbit collectively together with the unstable satellite 
by a uniform approximation discussed in ||. For k < £;( tan ) = 1.853, this approximation 
considers the central orbit as isolated and incorporates the satellites collectively j^, ||, |3| . It 
fails close to the bifurcations and only regains validity far away. There the orbits can also be 
treated individually and the stationary-phase approximation becomes justified. 

From the comparison of figs. ^| and |] it is apparent that the interval where all three orbits 
have to be treated collectively is quite large. It follows that unfoldings of bifurcations of higher 
co-dimension are relevant for semiclassical studies even if one varies only a single control 
parameter. This can be understood recalling that there exists a semiclassical measure of 
vicinity of orbits, the difference of their actions in units of Planck's constant, AS/ {2-kK). This 
measure is unknown to classical mechanics such that there bifurcations at different parameter 
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values do not feel each other so as to necessitate a qualitatively different treatment. 

In conclusion, in this Letter we studied configurations of periodic orbits that participate 
in sequences of bifurcations. It was shown that in semiclassical studies involving periodic- 
orbit theory such configurations have to be treated collectively. This was illustrated for the 
sequence of a period tripling and the tangent bifurcation of the satellite. A uniform collective 
contribution was given and tested for the kicked top. We found excellent agreement between 
semiclassical and exact peak heights in the action spectrum. 

In the language of normal-form theory, the case studied here is the unfolding of a period- 
three bifurcation of co-dimension two. Normal-form theory predicts similar scenarios for every 
period-71 bifurcation, and indeed further examples are under investigation. 
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